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Abstract 


Air-breathing proton exchange membrane (PEM) fuel cells provide for fully or partially passive operation and have gained much interest in 
the past decade, as part of the efforts to reduce the system complexity. This paper presents a detailed physics-based numerical analysis of the 
transport and electrochemical phenomena involved in the operation of a stack consisting of an array of vertically oriented air-breathing fuel cells. 
A comprehensive two-dimensional, nonisothermal, multi-component numerical model with pressurized hydrogen supply at the anode and natural 
convection air supply at the cathode is developed and validated with experimental data. Systematic parametric studies are performed to investigate 
the effects of cell dimensions, inter-cell spacing and the gap between the array and the substrate on the performance of the stack. Temperature and 
species distributions and flow patterns are presented to elucidate the coupled multiphysics phenomena. The analysis is used to determine optimum 
stack designs based on constraints on desired performance and overall stack size. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Increasing energy requirements around the world have stimu- 
lated the interest to develop new reliable clean energy to decrease 
the dependence on fossil fuels and to reduce atmospheric pollu- 
tion. To this end, proton exchange membrane (PEM) fuel cells 
are a strong alternative as a portable power source in widespread 
applications including automotive, laptop computers, cellular 
phones and other electronic devices owing to its low operating 
temperature, quick start, light weight and high power density. 
Recently, air-breathing PEM fuel cells have gained much interest 
due to their compact size and fewer components. Since oxygen 
is drawn directly from the atmosphere by natural convection, 
humidification and pressurization subsystems at the cathode are 
no longer required, which forms the basis of the reduced system 
complexity. 

Effective design of air-breathing PEM fuel cells requires 
a thorough understanding of the interrelated heat and mass 
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transport phenomena and electrochemical reactions in the cell, 
which has motivated the studies in the literature. Earlier experi- 
mental and numerical investigations on air-breathing PEM fuel 
cells focused on the performance evaluation and design of a 
single fuel cell. Several materials and fuel cell designs were 
explored [1-3], while analytical model and numerical simu- 
lations from physics-based models were used to predict the 
performance of the cell and to obtain detailed species and 
temperature distributions inside the cell [4,5]. A new design 
of an air-breathing fuel cell was reported by the authors [6], 
who also presented detailed computational modeling of the 
cell performance. The design consists of a dual fuel cell 
cartridge with a common hydrogen flow channel and two 
exposed cathode surfaces in communication with the ambient 
air. 

The studies on individual air-breathing fuel cells have shown 
the practical viability of the designs. Meeting the power require- 
ments of practical applications requires configuring multiple 
individual air-breathing cells connected either in series or in 
parallel as needed to form a stack. To the best of the authors’ 
knowledge, relatively scant work on air-breathing fuel cell stacks 
has been reported in the open literature. Chu et al. [7,8] reported 
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Nomenclature 


Ef 


Sv 


A 
e 
Q 
~ 


a~n 


water vapor activity 

fuel cell active area (m°) 

fixed charged site concentration in membrane 
(mol m~?) 

specific heat at constant pressure (J (kg Ky 
diffusivity of species i (m? s7!) 

cell voltage (V) 

stack voltage (V) 

Faraday’s Constant (Coulomb mol`!) 

gravity vector (m? s7!) 

vertical gap (m) 

stack height (m) 

average current density (A m7?) 

transfer current density (A m73) 

permeability of porous media or equivalent 
porous media (m?) 

electrokinetic permeability of membrane (m?) 
hydraulic permeability of membrane (m7) 
thermal conductivity (W (m K)~ 1) 

stack length (m) 

molecular weight of species i (kg s7!) 

flow rate (lpm) 

pressure (Pa) 

saturation pressure of water vapor in Eq. (11)(bar) 
relative humidity 

source terms of governing equations in Table 1 
temperature (K) 

superficial velocity vector (ms~!) 

volume of the catalyst layer (m?) 

position coordinate across the cell thickness/the 
membrane 

mole fraction of species i 

mass fraction species 7 

position co-ordinate along the chamber 

charge of sulfonate site in the membrane 


Greek letters 


transfer coefficient 

porosity 

volume fraction of membrane in the catalyst layer 
phase potential (V) 

surface over potential (V) 

water content 

fluid density (kg m~?) 

electrical/protonic conductivity (1 (ohm m)~!) 
stress tensor (kg (m s)~2) 


Subscripts/superscripts 


anode 

average 
cathode 
effective value 
fluid 
hydrogen 


i species 

m membrane 

Oo oxygen 

ref reference 

s solid 

t thermal 

Ww water 

x,y component values along the x and y coordinate 
directions 


experimental tests combined with numerical analysis based on 
empirical correlations to evaluate the steady-state performance 
of an air-breathing PEM fuel cell stack under different oper- 
ating conditions, while Morner and Klein [9] presented the 
dynamic behavior of an air-breathing fuel cell stack and the 
influence of temperature, humidity and air flow on the stack 
performance. Most recently, Litster and Djilali [10] imple- 
mented an empirical correlation for membrane conductivity 
to improve the accuracy of the modeling for dry membrane 
condition and the model was well validated with experimen- 
tal data. A comprehensive parametric study was conducted to 
investigate the effects of design specifications and operating 
conditions on the cell performance. The fundamental differ- 
ence between the conventional and the air-breathing fuel cell 
stacks is that in a conventional PEM fuel cell stack, the stack 
geometry and dimensions affect only the heat transfer charac- 
teristics, while in an air-breathing stack the internal and external 
stack designs affect both the heat transfer and the gas supply 
directly. 

A comprehensive analysis of the governing multiphysics phe- 
nomena is imperative for effective stack designs and forms the 
focus of the present work. Considering the dual-fuel-cell car- 
tridge design, described in Ref. [11], the study focuses on the 
analysis of a stack of such cartridges. The stack is configured 
as a vertical array of equally-spaced cartridges on a horizon- 
tal substrate, with a clearance between the substrate and the 
bottom face of the cartridge. The inter-cartridge spacing, the 
gap between the cartridge and the substrate, as well as the car- 
tridge dimensions affect the natural convective flow patterns at 
the cathode faces that influence the stack performance. This 
paper presents numerical studies with the objective of system- 
atically analyzing the effects of the stack height, the gap above 
the horizontal substrate, and spacing between cartridges on the 
performance of an air-breathing PEM fuel cell stack. A com- 
prehensive two-dimensional physics-based model is developed 
based on the model for a single cell presented in Ref. [11], and 
extended to include the effects of adjacent cartridges as well as 
the substrate underneath the cartridges. The numerical model is 
first validated with experimental data, which forms the basis for 
identification of the optimum stack parameters for maximizing 
performance subject to constraints on overall stack size. The 
numerical model adopted in the present study is presented in 
the next section, and the results of the studies are discussed in 
Section 3. 
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Table 1 
Source terms in the governing equations for the various regions of the fuel cell 
Region Se Sv ST Si Ss, Sm 
Gas channels 0 0 0 0 0 
Diffusion layers 0 — E eeft V 0 0 0 
and metal 
meshes 
Catalyst layers Anode: — Mp, - teen + ez Vebm Anode: jana + i Anode: — J Mm, (H2), 0 Anode: Ss = — ja < 0, 
" | (02), 0 (H20) Sm = +ja > 0 
Cathode: f Cathode: — jene + am Cathode: 0 (H2), | Cathode: S; = +jc > 0, 
=j Mo, + 3 Mmo " — {p Mo, (O2), $ Mmo Sm = —Jje < 0 
; (H20) 
Membrane 0 = EtV + erce F Vom “ar 0 0 


Active Natural Convection 
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Fig. 1. Schematic of an air-breathing PEM fuel cell stack. 
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2. Mathematical model 


Fig. 1 shows a schematic description of an air-breathing PEM 
fuel cell stack, consisting of an equispaced array of air-breathing 
fuel cell cartridges oriented vertically on a substrate. The stack 
parameters considered in this study are shown in Fig. 1: stack 
length, Ls; stack height, H; spacing between two adjacent fuel 
cell cartridges, 2b; and vertical gap between the stack and the 
bottom substrate, h. 


> 


uonoonu0D, 
jenjen 


Cartridge 


Cartridge 


Each air-breathing fuel cell cartridge in the stack consists 
of two single cells sharing a common hydrogen flow chamber 
on their anode sides, while the cathode sides of the two cells 
are exposed to ambient air. Fig. 2 shows an exploded schematic 
illustration of a two-dimensional cross section of the assembled 
cartridge spanning across the thickness of the cell (x-direction) 
and along the length of the hydrogen chamber (y-direction). Each 
cell in the cartridge comprises a layered structure of seven com- 
ponents as identified in one of the cells in Fig. 2: (1) an outermost 
frame with ribs to provide overall support; (2) a cathode metal 
mesh which acts as cathode current collector; (3) a cathode 
gas diffusion layer (GDL); (4) a membrane electrode assem- 
bly (MEA), consisting of a proton transport membrane and two 
catalyst layers on its surfaces where electrochemical reactions 
take place; (5) an anode gas diffusion layer; (6) an anode metal 
mesh which acts as anode current collector; (7) a common anode 
hydrogen flow chamber shared by the two cells. The component 
layers of the two cells are bolted together to form a cartridge. 
The two individual cells can be connected in parallel or series to 
meet the design requirements. The frames and the metal meshes 
serve to increase the stiffness of the cartridge and minimize the 
deformation of the fuel cell, and correspondingly decrease the 
interfacial contact resistances inside the cell. Due to the sym- 
metry along the mid-line of the hydrogen channel and that of 
cartridge spacing, the computational modeling is based on the 
half domain of the cartridge, identified by the thick dashed line 
in Fig. 2. The bottom line of the domain is determined by the 


1. Frame 
2. Anode metal I 
mesh 
3. Anode GO" 
4. MEA j 
5. Cathode GDL Cartridge 
6. Cathode metal 
mesh 
7. H, chamber 


Fig. 2. Schematic illustration of a two-dimensional section of an air-breathing PEM fuel cell cartridge in the stack. 
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edge of the substrate and the top line is set to ensure the enough 
extension region to simulate the natural convection flow outside 
the cathode. 

During the stack operation, hydrogen gas (H2) is supplied 
into the anode chamber of individual cartridge at a constant 
flow rate through a manifold, while the cathode surfaces are 
exposed to the atmosphere and the required oxygen (O2) is 
drawn into the cell from the ambient air around the stack by 
natural convection and diffusion. Therefore, the air flow pat- 
tern around the cathode side significantly affects the air supply 
at the cathode, which is determined by the geometry parame- 
ters of the stack including the cartridge spacing, 2b, the vertical 
gap, h, and the stack height, H. More fresh air will be provided 
to the open cathode with increase of the cartridge spacing and 
the vertical gap, which improves the air transport into the cell 
and correspondingly enhances the cell performance at a cost of 
larger size of the stack. Increasing the stack height increases the 
flow path of air; however, the air concentration will decrease 
along the height due to the reaction consumption, which may 
cause insufficient oxygen availability toward the top of the 
stack. 

To investigate the cell performance in the stack, a computa- 
tional modeling is developed to describe the physical phenomena 
inside a single cell in a stack. The gas transport through the 
domain is governed by the coupled continuity, momentum, and 
species conservation equations; the water transport through the 
membrane from anode to cathode is modeled as source terms 
representing the osmotic drag, and is included in the momen- 
tum equations for the catalyst layers and the membrane; and the 
electron transport through the external circuit is governed by 
solid phase potential equation in the catalyst, gas diffusion lay- 
ers and metal mesh at anode and cathode sides, while the proton 
transport in the MEA is expressed as membrane phase potential 
equation, in which the electrochemical reactions contribute as 
the source terms based on the Butler-Volmer expressions. The 
heat released from the electrochemical reactions and Joule heat- 
ing heats up the stack and the surrounding ambient air and are 
included in the energy equation to solve for the temperature dis- 
tribution inside the cell and between the cartridges, which drives 
the natural convection flow around the cell. 

The assumptions for the present model follow those rou- 
tinely adopted in the previous studies, namely: (1) the species 
are treated to be ideal gas mixtures; (2) the flow is assumed 
as laminar in the entire model domain; (3) liquid water is 
not considered due to the relatively low current density char- 
acteristic of air-breathing cells in comparison to conventional 
PEM fuel cells with forced cathodic flow; (4) the gas diffu- 
sion layers, catalyst layers and membrane are considered to 
be isotropic and homogeneous; (5) the membrane is consid- 
ered impermeable to gas flow and crossover phenomena are 
not considered; (6) contact resistances are not considered; (7) 
Joule heating is only considered in the membrane and catalyst 
layers. 

The coupled system of governing equations presented here 
is formulated to describe the entire half-domain in Fig. 2, with 
the phenomena specific to each region expressed in terms of 
appropriate source terms. This provides a unified formulation 


as summarized below with the specific source terms listed in 
Table 2: 


V-(pV) = Se (1) 
V-(oVV) =—Vp+V-t+ pet Sv (2) 
V- (pVyi) = V- (PDV yi) + Si (3) 
V - (pcpVT) = V - (KeEVT) + Sr (4) 
v- (ove) +S;=0; j=sm (5) 


where p is the fluid mixture density; V the superficial (volume- 
averaged) velocity vector with components u in the x-direction 
and v in the y-direction to ensure mass flux continuity at the 
interface between porous and non-porous media of the domain; 
p the pressure; T the stress tensor; g the gravity vector, which is 
Oms~? in the horizontal direction and 9.8 m s7? in the vertical 
direction; y; represents the i-th species mass fraction; cp is the 
fluid specific heat at constant pressure; Tis the temperature; and 
és and m denote the solid phase and membrane phase potentials, 
respectively. 

The source terms in the governing equations represent phys- 
ical phenomena specific to the various regions as follows: S¢ is 
mass consumption or generation terms in the two catalyst layers; 
Sy represents the pressure gradient due to flow through porous 
media, and includes both pressure gradient due to flow through 
porous media and electroosmotic force in the MEA since pro- 
tons are driven by osmotic force and transported from the anode 
catalyst layer to the cathode catalyst layer with water; S; stands 
for the consumption of H2 in the anode catalyst layer, and con- 
sumption of O2 and the generation of H20 in the cathode catalyst 
layer; Sr represents the heat generated in the catalyst layers due 
to reactions and Joule heating in the MEA. Joule heating in the 
solid phase is neglected since the electrical conductivity of the 
solid phase is much larger compared with the ionic conductiv- 
ity in the MEA; the source term S; denotes current generation 
owing to the electrochemical reactions in the two catalyst lay- 
ers, expressed by transfer current densities. The expressions for 
these source terms are summarized in Table 1. 

The species conservation equation (Eq. (3)) is applied to solve 
for the mass fractions of only n — 1 species components, and the 
last component is calculated by: 


n-1 


yn =1-S yi (6) 


i=l 


In the present study, since pure H2 is supplied into the car- 
tridges of the stack, the mass fraction of H2, O2 and H20 — yp,, 
yo, and yp,o — are obtained from Eq. (3), and that of N2, yn,, 
is calculated using Eq. (6). 

The governing equation set involves several parameters that 
are determined using appropriate models as summarized below: 


(1) The effective diffusivity of individual species in each region, 
DË, is modeled using Bruggeman correlation based on the 
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Table 2 
Physical and geometry parameters used in the numerical simulations 
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Parameter Symbol 


Length of Hz chamber, [cm] H 
Metal mesh (screen) width, [mm] 
Gas diffusion width, [mm] 
Catalyst layer width, [mm] 
Membrane width, [mm] 
Hz chamber width, [mm] 
Permeability of diffusion layer [m7] k 
Gas diffusion porosity £ 


Anode reference exchange current density 
[A/m>°] 


Cathode reference exchange current 
density [A/m*] 


Hydraulic permeability of membrane [m ? ] kp 

Electrokinetic permeability [m°] ky 

H20 viscosity in the membrane [kg/(m-s)] 
Anode/cathode transfer coefficient 

Fixed site charge Zt 

Membrane water porosity Em 

Volume fraction membrane in catalyst layer Eme 

Fixed charge concentration [kmol/m* ] Cr 


H, diffusivity in catalyst layer and membrane 
[m*/s] 


H) diffusivity in gas [m/s] 
O; diffusivity in gas [m7/s] 
H,O diffusivity in gas [m7/s] 


Oz diffusivity in catalyst layer and 
membrane [m/s] 


Value 


1.68/5.00 
0.5 
0.3 
0.016 
0.018 


2.97 x 10” 
0.776 


1.0x 10° 


3000 


1.58 x 1078 
1.13x 10° 
8.91 x 107 


0.5/2.0 


-10 


2.59 x 10° 


11x 107 
3.2x 10° 
7.35x 10° 


8.328 x 107 


Source [Ref.] 


Experiment [6] 
Experiment [6] 


J 
| 
5] 


[17] 


[14] 


effects of porosity and tortuosity in the porous media [12]: 


1.5 
T Pref 
i eff “i,ref i p 


(7) 


in which, €eg¢ is the porosity in the gas diffusion layers, 


(2) 


(3) 


catalyst layers and metal mesh and the equivalent porosity 
in the membrane, and Dj rer is the diffusivity of the species 
i at the reference temperature and pressure: Tref = 353 K 
and Pret = 10, 1325 Pa, listed in Table 2. Note that in the 
non-porous region, De defaults to 0 in a solid (€e¢¢ = 0) or 
the temperature- and pressure-dependent individual species 
diffusivity in an open channel (€e¢¢ = 1). 

The effective thermal conductivity, kf, is expressed as fol- 
lows in the porous regions [18]: 


kef = eeppki£ + (1 — Eetf)kr,s (8) 


where k; f and krs are thermal conductivities of fluid and 
solid, respectively. 

The effective electrical conductivity of solid phase, ot, is 
expressed as a function of the effective porosity, £eff, and the 
electrical conductivity of the solid material, o,, as a” = 


(1 — £eff)os; the effective membrane ionic conductivity in 
the catalyst layers can be expressed as: 


eff 1.5 
Om = En Om = ehi 


5(0.51392 — 0.326) 


1 1 
x exp 1268 (= = 7) 


where £mc is membrane volume fraction in the catalyst layer, 
and the ionic conductivity in the membrane phase, om, is 
expressed as the function of temperature, T, and the hydra- 
tion parameter, A, defined as water molecules per sulfonate 
group, in the membrane [13]. The hydration parameter, A, is 
estimated by the following empirical expression [13], such 
that a value of A larger than 14 indicates the presence of 
over-saturated water vapor in the membrane: 


(9) 


_ J 0.043 + 17.81a — 39.85a* + 36.0a°, O<a<1 
~ | 14414 - 1), l<a<3 
(10) 


where a is water vapor activity, expressed as a = Xw p/ Psat» 
in which, xw is water mole fraction, and Psat is the 
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temperature-dependent saturation pressure of water [13]: 
log Psat = —2.1794 + 0.02953(T — 273.15) 

—9.1837 x 107°(T — 273.15) 

+1.4454 x 1077(T — 273.15) (11) 
In the above expression, the saturation pressure, Psat, is in 


bar, and the temperature, T, is in Kelvin. 


As indicated in Table 1, the source term Sj, j= s,m, in 
the potential equations can be determined by transfer current 
densities, ja and je, governed by the Butler—Volmer equations: 


. ef *Ho a ep a, F i i a, F à 
A= Cee RT” RT” 
(12 
. zef X% e ack e acF 
Je = Je XOp,tef P RT Nc p RT Nc 


(13) 


in which, the subscripts a and c denote anode and cathode, 
respectively, jref is the reference volumetric transfer current den- 
sity, æ the transfer coefficient, xy, and xo, are the mole fractions 
of Hz and O3 in the catalyst layers, respectively, F is the Faraday 
constant, and 7 is the surface over potential defined as [16]: 


n = $s — bm — Ever (14) 


where Eyer is O at the anode and equals the temperature- 
dependent open cell voltage, Eoc = 0.00257 + 0.2329, on the 
cathode side, where the temperature, 7, is in Kelvin [14]. 

The boundary conditions associated with the governing equa- 
tions set are shown in Fig. 3. The mid-line of the hydrogen 
chamber and the mid spacing between cartridges are both 
assumed to be lines of symmetry due to geometry symmetry 


P, =ConSt. Dp, = Pa T; =T Vig = Yar Pe = Px 
y, ə- EL Du 
CEE 0 7 =0 
ey x | | & 
T ohn > OT. 
9 ae | <=0 
ey D DE ee E 
<3 5 SS | 
| ee 3 mae s 
| g% = o | 
oe ie 
| ee FS | 
KS x 
=: : | 
symmet ric H. 2 oe s Air | symmet ric 
¢, =0 vi oe ¢, = const. | 
| Se 
Kee 
PE, | se N 
| ee | ‘Model 
v = const. $0909 š 
ass | Domain 
Yy, = Const. Set | 
p Ô ð 
Yo, =0 | On =0 ge 
ox | ox 


Yy o = CONSI. j= -—---<-- 


T, =const. Pe = Pas T, =T 5 Vie = air heated wall (q" „= const.) 


Fig. 3. Schematic of the model domain and the associated boundary conditions. 


and nearly identical operation of each cartridge in the stack. In 
the present study, the substrate is considered to be a DC-DC 
power converter generating constant heat; therefore, a constant 
heat flux of 0.025 W m~*was specified at the bottom surface of 
the domain. The interface conditions at any two adjacent layers 
in the domain are implicitly satisfied since the governing equa- 
tions are formulated and implemented for the entire region as 
a single domain. Accordingly, the boundary conditions for the 
solid phase potential equation, namely, ¢, = 0, at the interface 
of the hydrogen chamber and the metal mesh and ¢, = Ee at 
the cathode surface are implemented at the outer left and right 
borders of the domain with infinite electrical conductivity in the 
hydrogen chamber region and in the extended ambient region. 
Since there is no proton transport outside the MEA, the boundary 
conditions for the membrane phase potential are zero proton flux 
at the interfaces of the catalyst layers and gas diffusion layers, 
which is implemented as zero ionic conductivity in the regions 
except the MEA. The boundary conditions at hydrogen inlet 
are specified velocity components, u, v, specified mass fractions 
corresponding to humidified hydrogen, yy, a, YH20,a; Specified 
anode temperature, T,; while those at the hydrogen outlet are 
as follows: specified anode pressure, pa, and zero gradients for 
other flow quantities. 

The model described above was implemented in a commer- 
cial finite-volume based computational fluid dynamics solver, 
Fluent, using the physical and geometry parameters listed in 
Table 2. Several user-written subroutines were embedded to 
describe mass consumption/generation, heat generation, water 
transport, electrochemical reactions and the potential equations 
as well as the effective properties for specific regions. A segre- 
gated implicit method was used to solve the coupled governing 
equations and the semi-implicit method for pressure-linked 
equations consistent (SIMPLEC) algorithm [19] was used for 
pressure correction to ensure mass conservation. No commer- 
cially available fuel cell module from Fluent was used. In the 
present study, the cells in the stack are considered to be con- 
nected in series and the performance of each cell is assumed 
identical in the stack. Therefore, the individual cell voltage, 
Ecen, which is one of the boundary conditions of the numeri- 
cal modeling, is calculated as the stack voltage, Estack, divided 
by the the number of the cells in the stack. Since each cartridge 
consists of two cells, the cell voltage, Ee, is equivalently the 
stack voltage divided by twice the number of cartridges in the 
stack. According to the boundary conditions shown in Fig. 3, 
the coupled set of governing equations was solved iteratively 
to obtain the velocity, pressure, temperature, species and two 
potential distributions, from which the volumetric transfer cur- 
rent densities, ja and je were calculated through Eqs. (12) and 
(13). The average current density of the fuel cell, Z, is calculated 
by integrating the volumetric transfer current densities over the 
entire volume of the catalyst layer, and then dividing by the 
active area in the y—z plane (shown in Fig. 1), A, as: 


1 1 
I= 5 | av = = | av (15) 


Therefore, the polarization curve is constructed by solving for 
the average current density, Z, for different values of cell volt- 
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Stack Voltage, E tack [V] 


Experiment [6] 
Numerical Simulation 
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2,500 


3,000 


Fig. 4. Validation of the numerical simulations for the stack with experimental 
polarization curve data from Ref. [6]. 


age, Ecel. The convergence criterion is set as that the maximum 
scaled residual of all the variables being solved is less than 10~°. 
The numerical simulations were used to explore the cell perfor- 
mance over a range of geometry parameters and the results are 
presented in the following section. 


3. Results and discussion 


The two-dimensional physics-based model presented in the 
previous section is first validated with experimental data from an 
air-breathing PEM fuel cell stack [6] consisting of cartridges as 
described in Fig. 2. Also, the numerical model was validated with 
the experimental data on the single cell performance for different 
operating conditions in Ref. [11]. Fig. 4 presents the comparison 
of polarization curves from the numerical simulation denoted as 
the solid line, and the experimental data from Ref. [6] described 
by the solid markers, for a stack consisting of five fuel cell car- 
tridges, each of height, H = 5.0cm, with cartridge spacing of 
2b = 20.0 mm, and vertical gap of h = 15.0mm. The ten cells 
in the five cartridges are connected in series, and the current is 
collected through a square active area of 25 cm”. The stack oper- 
ating conditions are as follows: room temperature at the anode, 


WMMMM«Mab 
10.0 mm 


h=0.0 mm 0.5 mm 2.0mm 


Ta = 298 K, anode pressure, pa = 1.0 atm, zero anode relative 
humidity, RH, = 0.0, and constant anode flow rate for each car- 
tridge, ma = 0.15 lpm, while the cathode boundary conditions 
are as shown in Fig. 3. For the following discussion in this sec- 
tion, the above operating conditions are referred to as the base 
operating conditions. The stack voltage is calculated by multi- 
plying the cell voltage from the numerical simulation by ten, 
the number of cells in the stack used in the experimental study 
[6]. It is seen in Fig. 4 that the stack voltage, Estack, decreases 
monotonically with increase of the average current density, Z, as 
expected due to ohmic loss in the membrane, and activation and 
concentration losses in the catalyst regions. The model predic- 
tions in Fig. 4 show a close agreement with experimental data, 
which forms a reliable basis for the use of the present numerical 
model to conduct systematic parametric studies. Since the stack 
voltage is directly related to each individual cell voltage by the 
number of cells, in the following discussion, the performance of 
a single cell in the stack is investigated. 

Fig. 5 presents the temperature contours and correspond- 
ing flow patterns in the modeling domain for different vertical 
gaps and cartridge spacings, and for a fixed cell voltage of 
Ecen = 0.5 V and stack height of H = 5.0 cm under the base 
operating conditions. The temperature variation is represented 
by the dashed contour lines with corresponding temperature 
values in Kelvin and the stream lines are represented by the 
solid lines with arrows indicating the flow direction. The effects 
of vertical gap on temperature distributions and flow patterns 
are shown in Fig. 5a, where the cartridge spacing is fixed at 
2b = 20.0 mm. It is observed that the maximum cell tempera- 
ture is approximately 303 K and the temperature variation inside 
the fuel cell is relatively small, around 1 K. The average tem- 
perature in the inter-cartridge space is seen to decrease with 
increasing vertical gap, which can further be explained from air 
flow patterns for different vertical gaps. If the bottom of the stack 
touches the top of the substrate, i.e., h = 0.0, air flow forms a 
recirculating pattern in the half domain of the inter-cartridge 
space (shown in the leftmost frame in Fig. 5a), and therefore, 
only a limited amount of fresh ambient air could enter the circu- 
lation inside the inter-cartridge space from the top of the stack. 
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Fig. 5. Effects of (a) vertical gap, h, and (b) cartridge spacing, 2b, on the temperature distributions and flow patterns of an individual cell in the stack at H = 5.0 cm. 
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Fig. 6. Effects of (a) vertical gap, h, and (b) cartridge spacing, 2b, on the temperature distributions and flow patterns of an individual cell in the stack at H = 1.68 cm. 


As vertical gap increases, ambient air flows from the vertical 
gap and ascends along the entire height of the cartridge to pro- 
vide natural convective cooling which leads to lower average 
temperature in the cartridge space. If the vertical gap is larger 
than a certain critical value, the temperature in the middle of car- 
tridge space is reduced to room temperature (lower than 299 K 
in the plot), signifying that the convective flow pattern for an 
individual cell in the stack is no longer affected by the adjacent 
cells. 

The effects of cartridge spacing on temperature distributions 
and flow patterns are shown in Fig. 5b, where the vertical gap 
is fixed at the largest value studied, h = 10.0 mm, and the spac- 
ing is varied from 2b = 4.0mm to 2b = 12.0mm. It is seen 
that the temperature distributions and flow patterns are similar 
for the different cartridge spacings. However, the temperature 
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value is seen to be higher for small cartridge spacings due to 
the proximity effects of the adjacent cells. Fig. 6a and b show 
the temperature distributions and flow patterns under different 
vertical gaps and cartridge spacings for a shorter stack height of 
H = 1.68cm. The temperature distributions and flow patterns 
have the same qualitative trends as those for the larger stack 
height of H = 5.0cm in Fig. 5. The results indicate that varia- 
tion in the stack height has the least effect on the thermal and 
fluid transport characteristics of the fuel cell stack, as opposed to 
the cartridge spacing, 2b and vertical gap, h, which significantly 
affect the temperature distribution and the flow patterns in the 
fuel cell stack. 

The effects of vertical gap, h, and cartridge spacing, 2b, on the 
cell performance for the stack height, H = 5.0 cm, are shown in 
Fig. 7. Fig. 7a presents polarization curves of an individual cell in 
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Fig. 7. The performance of an individual cell in the stack under different vertical gaps, h, and cartridge spacings, 2b, at H = 5.0 cm: (a) the polarization curves for 
different h; (b) the current density variation with h for different cell voltages; (c) the polarization curves for different 2b; (d) the current density variation with 2b for 


different cell voltages. 
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the stack with different vertical gaps, h, at a fixed cartridge spac- 
ing of 2b = 20.0 mm. Note that with increase of h from 0.0 mm 
to 0.5 mm, the cell performance is significantly improved, while 
with increase of h from 0.5mm to 2.0mm, the cell perfor- 
mance only slightly improves, and with further increase of h 
from 2.0mm to 10.0mm, the performance remains the same, 
which indicates that beyond a certain clearance between the 
cartridge and the bottom substrate, the cell performance cannot 
be enhanced further with increase of h. This trend is corrobo- 
rated by the flow patterns in Fig. 5a, which shows that increasing 
h to a non-zero value completely changes the flow patterns in 
the inter-cartridge space, whereby fresh ambient air is continu- 
ously drawn through the inter-cartridge space from the bottom 
to the top and results in enhanced oxygen transport into the 
cell for reactions. This pattern is improved with increasing of h 
from 0.5 mm to 2.0 mm in Fig. 5a; however, further increase of 
h results in no significant difference in the flow patterns. Cor- 
respondingly, the cell performance is seen to be only slightly 
improved for h = 2.0 mm to h = 10.0 mm since the cell reaches 
oxygen transportation limitation and no more improvement in 
the cell performance can be obtained from increasing the oxygen 
availability at the cathode. Fig. 7b shows the influence of h on 
the average current density for different values of Ecen: 0.2 V, 
0.3 V, 0.50 V, 0.70 V and 0.85 V. It is seen that for the larger cell 
voltages of Ecen = 0.85 V and 0.70 V, the average current den- 
sity is nearly invariant with the vertical gap, h, due to the small 
oxygen consumption at these cell voltages and the relatively 
sufficient oxygen supply even for a small clearance or no clear- 
ance. With decreasing cell voltage to Ece = 0.50 V, 0.30 V and 
0.20 V, the current densities increase sharply with increase of h, 
and saturate for h exceeding a minimum value, Amin. As shown 
in Fig. 5a, increasing h from 0.0 mm improves the air supply at 
the cathode: for h < hmin, the fuel cell performance is limited 
only by the air supply and improving the air supply results in 
enhanced current density; however, further increase in h does 
not affect the current density owing to the air supply exceeding 
the need corresponding to the cell voltage. Furthermore, it is 
noted in Fig. 7b that the minimum height, hmin, above which the 
current density saturates increases with decreasing cell voltage. 

The polarization curve of a single cell in the stack for different 
cartridge spacing, 2b, and at a fixed vertical gap, h = 10.0 mm, 
is presented in Fig. 7c. The polarization curves overlap at the 
higher voltages or equivalently, the smaller current densities, 
indicating the negligible influence of cartridge spacing since 
the air supplied through the smaller inter-cartridge space is still 
sufficient for stack operation at the smaller current densities. 
As the cell voltage decreases, the current density decreases 
with the inter-cartridge spacing, and the effect is seen to be 
especially pronounced for 2b = 4.0 mm. A small inter-cartridge 
space reduces the air flow rate from the bottom to the top in the 
space, as shown in Fig. 5b, and results in insufficient availabil- 
ity of oxygen for the electrochemical reactions, correspondingly 
impairing the cell performance. Fig. 7d depicts the influence of 
inter-cartridge spacing on the current density for different cell 
voltages. It is observed that current density initially increases 
with increase of cartridge spacing at different cell voltages since 
oxygen supply is increased in proportion to the inter-cartridge 


spacing. However, as the cell reaches its oxygen transport limit 
corresponding to each cell voltage, no additional benefit is real- 
ized by increasing the cartridge spacing and the current density 
is seen to saturate beyond a certain critical spacing. The criti- 
cal cartridge spacing is noted in Fig. 7d to increase from about 
1.2mm for Ecen = 0.85 V to 8.0mm for Ecen = 0.20 V. 

Fig. 8 presents the effects of the vertical gap and the inter- 
cartridge spacing on the cell performance in the stack of a shorter 
height of H = 1.68 cm. The polarization curves for different ver- 
tical gap, h, and the fixed inter-cartridge spacing, 2b = 20.0 mm, 
are presented in Fig. 8a, where it is shown that increasing the ver- 
tical gap from 0.0 mm to 2.0 mm improves the cell performance 
at low cell voltages since the flow patterns are improved and 
more air is supplied from the clearance as shown in Fig. 6a. Fur- 
ther increase in the vertical gap does not affect the performance 
due to the sufficient air supply relative to the transport limitation 
inside the cell. The same trend is also reflected in Fig. 8b, where 
for Ecen = 0.35 V and 0.20 V, the current densities are improved 
significantly by increasing h in the range of 0.0-2.0 mm. Com- 
paring Figs. 7b and 8b, it is seen that for low cell voltages in 
the range of 0.5-0.2 V, the increase in the current density with 
h for H = 5.0cm is larger than that for H = 1.68cm, which 
indicates that the performance of the large cell is more sensitive 
to the change of the clearance, h. The effects of inter-cartridge 
spacing, 2b, on the cell performance in the stack are presented 
in Fig. 8c and d, where it is observed that as in Fig. 7c and 
d, there is a critical inter-cartridge spacing beyond which the 
current density is relatively invariant with the spacing. Fig. 8d 
further shows that the critical cartridge spacing increases with 
decreasing cell voltage, which corresponds to increasing oxy- 
gen consumption. Furthermore, a comparison of Figs. 7d and 8d 
reveals that the critical value of inter-cartridge spacing is larger 
for the larger cell size for a given cell voltage. Figs. 7 and 8 
show that a minimum (critical) vertical gap and inter-cartridge 
spacing are needed for maximizing stack performance, and that 
these values are a function of the stack height, H. 

To further illustrate the effects of the vertical gap and car- 
tridge spacing on the transport and electrochemical phenomena 
inside the fuel cell, it is worth examining the species distribu- 
tions on the cathode side. Fig. 9a and b present oxygen and 
water distributions, respectively, on the surface of the cathode 
GDL as a function of the vertical location along the stack height, 
y, for different vertical gaps, h, and H = 5.0 cm. In Fig. 9a, it 
is observed that with no vertical gap between the fuel cell car- 
tridges and the bottom substrate, i.e., h = 0.0mm, the oxygen 
mass fraction in the mixture is very low at about 0.04 along 
most of the surface of the cathode GDL, and increases rapidly 
from 0.04 to 0.20 at the upper end of the stack, reaching toward 
the oxygen mass fraction in the atmosphere of 0.23. This trend 
is consistent with air flow patterns in the inter-cartridge space 
shown in Fig. 5a and the corresponding performance in Fig. 7a. 
Without a vertical clearance, the only source for oxygen is from 
the top of the stack, resulting in very low oxygen mass fraction 
available for the electrochemical reactions. As the vertical clear- 
ance increases, a continuous flow path from the bottom to the 
top of the stack is established, which provides increasing oxy- 
gen supply to the cathode GDL surface. This trend is evident in 
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Fig. 8. The performance of an individual cell in the stack under different vertical gaps, h, and cartridge spacings, 2b, at H = 1.68 cm: (a) the polarization curves for 
different h; (b) the current density variation with h for different cell voltages; (c) the polarization curves for different 2b; (d) the current density variation with 2b for 
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Fig. 9. Effects of vertical gaps, h, and cartridge spacings, 2b on species distributions along the cathode surface for H = 5.0 cm: (a) oxygen mass fraction distribution 
and (b) water mass fraction distribution for different h; (c) oxygen mass fraction distribution and (d) water mass fraction distribution for different 2b. 
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the increasing oxygen mass fraction along the stack height in 
Fig. 9a for increasing values of h. For cases with h > 0.0mm, 
as the location along the channel, y, increases from the bottom 
to the top, oxygen mass fraction initially decreases due to the 
rapid electrochemical reactions, and as y approaches the top of 
the stack, the oxygen mass fraction increases to the value close 
to that in the atmosphere. As A increases, the magnitude of the 
total oxygen mass fraction variation along the channel is seen 
to decrease due to improved oxygen supply. It is also observed 
that as vertical gap increases from 2.0 mm to 10.0 mm, oxygen 
mass fraction distribution increases only slightly, indicating a 
stabilization of the flow pattern and the oxygen supply relative 
to the need for the electrochemical reactions. 

Water mass fraction distributions on the cathode GDL sur- 
face are shown in Fig. 9b corresponding to the conditions of the 
cell discussed in Fig. 9a. Using the water mass fraction distri- 
bution and the pressure and temperature distributions (Fig. 5), 
the relative humidity can be readily calculated. To compare 
with the oxygen mass fraction profiles shown in Fig. 9a, the 
water mass fraction instead of the relative humidity is presented 
and discussed here. In contrast to the oxygen mass fraction, 
the water mass fraction for a stack with no vertical clearance 
(h = 0.0 mm) is very high at the bottom of the cathode surface. 
As the flow patterns in Fig. 5a indicate, species including water 
are not effectively transported out of the cell, which leads to 
high water mass fraction on the surface at the bottom of the 
stack. Furthermore, the accumulated water practically blocks 
oxygen transport and results in degradation of the fuel cell 
and stack performance. As a vertical clearance is introduced 
by varying h from 0.0mm to 2.0mm, the water mass frac- 
tion on the cathode GDL surface decreases dramatically since 
the resulting flow pattern improves water transport, correspond- 
ingly reducing the concentration variation along the stack height. 
The water distribution also does not change appreciably as the 
vertical gap increases from 2.0mm to 10.0mm, following a 
pattern similar to that in the oxygen mass fraction profiles in 
Fig. 9a. The results in Fig. 9a and b indicate that a clearance 
of 2.0mm or more may be optimum from a stack performance 
consideration. 

Fig. 9c and d present the effects of cartridge spacing on the 
oxygen and water species distributions on the cathode GDL 
surface for the fixed vertical clearance of h = 10.0 mm. As car- 
tridge spacing, 2b, increases from the smallest value of 4.0 mm 
to 8.0mm, the oxygen mass fraction values increase, while the 
magnitude of the total variation in the profiles along the cath- 
ode surface decreases owing to the improved oxygen supply. 
However, increasing cartridge spacing from 8.0 mm to 20.0 mm 
does not change the oxygen distribution as the fuel cell reaches 
the species transport limitations and does not benefit from the 
oversupply of air made available through the increased spacing. 
Water mass fraction distributions for different cartridge spacings 
are presented in Fig. 9d, where it is observed that the trend in the 
profiles is the opposite to that of the oxygen distributions. As car- 
tridge spacing increases, water mass fraction along the surface 
decreases corresponding to the increase in the oxygen profile 
in Fig. 9c, and furthermore the variation in the value along the 
stack height also decreases. The profiles in Fig. 9c and d suggest 


that a spacing of 8.0mm may be optimum from a stack design 
viewpoint. 

Fig. 10 presents the oxygen and water distributions along 
the cathode GDL surface corresponding to the stack height of 
H = 1.68cm. Oxygen mass fraction distributions on the sur- 
face for different vertical gaps are shown in Fig. 10a, where 
it is observed that the variations follow a similar trend to the 
results for the stack height of H = 5.0cm in Fig. 9a, discussed 
previously. Fig. 10b shows the corresponding water mass frac- 
tions for different vertical gaps. A comparison of the results in 
Figs. 9a and 10a and those in Figs. 9b and 10b indicate that for 
the same vertical gaps, the values of the oxygen mass fraction 
are greater and the values of the water mass fractions are sig- 
nificantly lower for the smaller stack height, implying that for 
identical vertical gap and cartridge spacing, oxygen and water 
transport is improved by decreasing the stack height. Equiva- 
lently, the results imply that the value of the critical vertical 
clearance increases with increase in the stack height, H. Oxygen 
and water mass fraction distributions for different cartridge spac- 
ings are shown in Fig. 10c and d. As cartridge spacing decreases, 
the oxygen mass fraction decreases while water mass fraction 
increases, leading to a reduced performance of the stack. How- 
ever, the magnitude of the decrease is smaller compared with 
that for the larger stack height, indicating that the minimum 
required cartridge spacing also deceases with a decrease in the 
stack height, H. 

The performance analysis of an individual air-breathing cell 
in a stack, presented in the results so far, may be used towards the 
design of a stack, as illustrated in the discussion in the remainder 
of this section. The considerations underlying a stack design are 
generally (a) a prescribed voltage to be delivered by the stack, 
(b) maximization of the power density available of the stack, and 
(c) minimization of the overall stack dimensions. It was shown 
that increasing the clearance, h, and the inter-cartridge spacing, 
2b, improves the air flow patterns surrounding a cartridge and, 
in turn, the performance of the cartridge. However, increasing 
h and 2b also lead to increasing the dimensions of the stack. 
Similarly, a given voltage requirement from a stack may be real- 
ized by including many closely-spaced cartridges, each with 
a lower cell voltage due to the constricted flow patterns in the 
narrow inter-cartridge spaces, or a few appropriately spaced car- 
tridges, each with a larger cell voltage enabled by the improved 
flow between the cartridges. In any case, a minimum number of 
cartridges is required to meet the voltage requirements. Thus, 
the stack design parameters to be determined from the analysis 
are: (1) the clearance of the cartridges from the substrate, h; (2) 
the (minimum) number of cartridges, N, required; (3) the inter- 
cartridge spacing, 2b. For the discussion below, the cartridge 
height, H, is considered to be 5.0 cm, although the methodology 
outlined is equally applicable to any other value of H. 

Based on the variation of the current density with the vertical 
clearance, h, for the different cell voltages in Fig. 7b, the clear- 
ance, h, is taken to be 2.0 mm from the observation that a larger 
clearance serves no useful effect on the performance improve- 
ment of the cell, whereas undesirably increases the overall stack 
height. The polarization curve corresponding to a clearance of 
2.0mm, such as that shown in Fig. 7c, is cast in terms of a cell 
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Fig. 10. Effects of vertical gaps, h, and cartridge spacings, 2b on species distributions along the cathode surface for H = 1.68 cm: (a) oxygen mass fraction distribution 
and (b) water mass fraction distribution for different h; (c) oxygen mass fraction distribution and (d) water mass fraction distribution for different 2b. 


voltage per unit stack length, E’, by dividing the cell voltage 
by a nominal substrate length, (b + 0.51), associated with each 
cell (half cartridge), where t is the cartridge thickness (fixed at 
6.5 mm in this study). Fig. 1 1a recasts the information contained 
in Fig. 7c in terms of a plot of the current density as function 
of the inter-cartridge spacing for four different values of the cell 
voltage per unit length, E’. The following observations are evi- 
dent from Fig. | 1a: all the lines converge to a current density of 
zero in the limit of the inter-cartridge space approaching zero, 
since physically the air supply at the cathode is cut off in the 
limit of zero spacing; for each value of E’, there are two values 
of cartridge spacing, 2b, that yield the same current density — 
one corresponding to operation at a higher cell voltage at a larger 
spacing (denoted as point A on one of the curves in Fig. 11a) 
and the other corresponding to operation at a lower cell voltage 
at a smaller cartridge spacing (denoted as point B in Fig. 11a) 
— although the total power density per unit stack length is iden- 
tical for the two designs; and for each value of E’, there is an 
optimum spacing that maximizes the current density output of 
the stack, and this optimum spacing as well as the associated 
maximum current density both increase with £’. 

Fig. 11b presents the power density variation with the inter- 
cartridge spacing, for the four different values of E’. The plots 
are obtained from Fig. lla by multiplying the cell voltage at 
every point in the lines. The results indicate an optimum spacing 
that maximizes the power density for each F’, and a comparison 
with Fig. lla reveals that the optimum spacing for maximum 
power density is larger than that maximizing the current density. 


This suggests that for a fixed voltage output from the stack, the 
design for maximum power density could be achieved with fewer 
cell cartridges (each operating at lower current density) than 
the number needed for maximizing the current density (which 
corresponds to a lower power density from the stack). Therefore, 
the design for maximum current density offers the largest power 
for a given stack length, while the design for the maximum 
power density provides for a lighter stack with fewer cartridges. 
The points of equi-current-density indicted as points A and B 
in Fig. 11a are also indicated in Fig. 11b to illustrate that the 
design with a larger inter-cartridge spacing (point A) has a higher 
power density compared to B with closely spaced cartridges, 
even though both designs offer the same current density and the 
same stack voltage for a given stack length. 

Figs. lla and b can be used as design charts for stack 
design as illustrated in the following example application which 
calls for a stack power, Pstack, and voltage, Estack, of 40 W 
and 8 V, respectively, and a maximum stack length, Ls, of 
80mm. The specifications correspond to E’ = Egtack/Ls = 
8/80 = 0.1 Vmm7' with a series connection of the cartridges in 
the stack. From the viewpoint of minimizing the stack weight, 
the design point is set as that for the maximum power density, 
denoted as point C in Fig. 11b. Based on Fig. 11b, for E’ = 
0.100 V mm™!, the maximum power density is 1700 W mat 
an inter-cartridge spacing of 2b = 3.5 mm. Therefore, the num- 
ber of cartridges needed in the stack is: N = Ls/(2b + t) = 
80/(3.5 + 6.5) = 8, which leads to the voltage for each cell, 
Even = Estack/2N = 8/(2 x 8) = 0.5 V. Also, for the same £’ 
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Fig. 11. Effects of cartridge spacing, 2b, on (a) current density and (b) power 
density distributions at H = 5.0 cm. 


and 2b, the current density for point C in Fig. 1la is deter- 
mined as J = 3450 A m2. Since the stack current is: istack = 
Pstack/ Estack = 40/8 = 5 A, the active area of the fuel cell unit 
is: A = igtack/I = 5/(3450 x 1074) = 14.5 cm’. Since the car- 
tridge height, H, is 5.0cm in Fig. 11, the corresponding width 
of the cartridge is: W = A/H = 14.5/5 ~ 3.0 cm. The vertical 
clearance, h, is taken to be 2.0mm since the stack performance 
is invariant above this value. The stack design is thus summa- 
rized as follows: L; = 8.0cm, H = 5.0cm, W = 3.0cm, h = 
2.0mm, 2b = 3.5mm and N = 8, and the operating conditions 
are Ecen = 0.5 V and I = 3450Am~. This design method- 
ology is applicable to any other specifications of application 
requirements and illustrates a physics-based approach to fuel 
cell stack design. 


4. Conclusions 


A comprehensive numerical analysis of a stack consisting of 
an array of air-breathing fuel cell cartridges was presented to 


illustrate the effects of the various stack and cell parameters, 
as well as the use of the analysis for stack design. It was shown 
that a critical minimum clearance between the fuel cell cartridges 
and the bottom substrate and a minimum inter-cartridge spac- 
ing were needed to enable the air flow patterns required for 
best performance of the stack. Furthermore, optimum spacings 
for two different objectives of maximizing current density and 
maximizing power density were developed and a comparative 
discussion of stack designs corresponding to these spacings was 
presented. The use of the analysis for design of a stack for target 
performance requirements and constraints was also illustrated. 
Overall, the study demonstrated the viability of air-breathing 
fuel cell for practical applications. 
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